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Abstract 

We compute the Casimir energy for a fermion chirally coupled to a prescribed pseudoscalar field 
in the form of the soliton of the sine-Gordon model. The classical equations of this system are not 
analytically solvable. We obtain the bound energies of the Fermi field using a numerical method 
called the relaxation method, and its continuum scattering wavefunctions by the use of the Runge- 
Kutta-Fehlberg of order 6. Then, we obtain a linear combination of the obtained scattering states, 
which is the simultaneous eigenstate of the Hamiltonian and parity. Using the obtained scattering 
wavefunctions and the parity eigenstates, we are able to calculate both the scattering phase shifts 
and the phase shifts of the parity eigenstates. The resulting phase shifts are consistent with both 
the weak and strong forms of the Levinson theorem. Having these phase shifts we can compute 
the Casimir energy of the system in two different ways. We can use the scattering phase shifts to 
obtain the Casimir energy, or equivalently add the Casimir energy for each parity, calculated by 
the phase shift for the corresponding parity. We plot the total Casimir energy and also the Casimir 
energy of each parity as a function of the parameters of the pseudoscalar field. In the graphs of 
the total Casimir energy and the Casimir energy of the positive parity as a function of the scale 
of variations of the pseudoscalar field, fi, there is a maximum occurring when a fermionic bound 
energy level with positive parity crosses the line of zero energy. Since no energy level with negative 
energy crosses the line E — 0, as we increase fi, there is no maximum in the graph of the Casimir 
energy of the negative parity. The limits fi — > and /i — > oo of the Casimir energy are the ones we 
expect. The graph of the Casimir energy is in general an increasing function of the value of the 
pseudoscalar field at spatial infinity, 8o- In the graphs of the total Casimir energy and the Casimir 
energy of each parity as a function of #o there are mild cusps when a bound state crosses E = 0. 
Moreover, considering a system consisting of a valence fermion in the ground state, we find that 
the most energetically favorable configuration is the one with a soliton of winding number one. 

1 Introduction 

In the original papers in 1948, Casimir [TJ [5] calculated the force between two uncharged perfectly 
conducting plates separated by vacuum and found that despite the lack of an external field, there 
exists a net force on the metallic plates and they attract each other. The Casimir effect arises from 
the change in the zero-point energy of the system. This change could happen when different boundary 
conditions are imposed or non-trivial spatial background fields such as solitons are present. 

The Casimir effect often is pointed out as an evidence of the reality of the quantum fluctuations of 
fields in the vacuum. Since the energy of the vacuum does not seem to be observable in any laboratory 
experiment, the appeal of the Casimir force as a macroscopic manifestation of the quantum zero-point 
fluctuations has drawn much attention over the decades that followed its discovery. Although there 
are a vast amount of works [3j SI [5] in which theorists refer to the Casimir effect as a good evidence 
of the reality of the vacuum fluctuations, recently some authors [5] believe that this effect gives no 
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support for the subject (reality of the vacuum energy). Their reason is that the Casimir force can be 
expressed entirely in terms of Feynman diagrams with external legs and without any reference to the 
vacuum energy. 

The Casimir energy and the resulting forces have been investigated for different fields in different 
geometries, as the boundary conditions of the system, [7]. Also, the majority of the investigations 
related to this effect concern the calculation of the Casimir forces between the boundaries in arbitrary 
spatial dimensions. This force could be either repulsive or attractive. The fact that the energy of a 
fluctuating field diverges both in the free space vacuum and in a vacuum in which a boundary condition 
or a background field is imposed, have been known for many years [5]. Many different regularization 
schemes or renormalization programs have been proposed to cancel these divergences when we subtract 
the energy of the free space vacuum from the energy of the vacuum in the presence of the boundary 
conditions or the background fields to find the Casimir energy [S]. 

One of the first experimental attempts to observe this phenomenon was conducted by Sparnaay 
[10] in 1958. Their measurement with parallel plates, though not in contradiction with the Casimir 
theory, could not be an evidence for the Casimir effect, due to large experimental errors. Quantitative 
experimental confirmation was just possible after the advent of the high precision probes in the 1990s. 
In 1997 LamoreauxQl] [T2] conducted the first successful measurement for the Casimir effect using a 
plate and a part of a sphere with a large radius rather than two parallel plates which would require 
accurate alignment to ensure they are parallel. Since then, many experimental investigations have 
measured the Casimir force in different shapes of the boundary conditions [13] . 

Recently, there have been many studies to use the Casimir force for practical applications. As device 
dimensions for Micro and Nano Electromechanical Systems (MEMS and NEMS) decrease, forces such 
as the Casimir forces which are normally neglected in macro systems have to be considered. The 
Casimir force could hinder MEMS and NEMS because it could cause the tiny elements in the device 
to stick together and helps us to design the wear-proof nano-scale systems [HI [Tj2 EH El [15] . 

As mentioned above, the zero-point energy can also be affected by the presence of non-trivial 
background fields which are usually chosen to be soliton. Also sometimes a very simple potential such 
as an electric potential well is chosen as the background field. This simple choice renders the problem 
of vacuum polarization and the Casimir energy exactly solvable and helps us to see more details |19j . 
Some authors compute the Casimir energy as the lowest order quantum correction to the mass of the 
soliton for models containing solitons such as supersymmetric solitons [2Ql [211 1221 [221 [2H [25] ■ For 
most of the models with solitons, the problem is not exactly solvable and the Casimir energy cannot 
be calculated directly. However, it can be calculated by the use of the indirect methods such as the 
phase shift method which relates the derivative of the phase shift with respect to the momentum to 
the spectral deficiency in the continuum states [501 [221 1221 12H EH ■ 

In the present paper we choose a Lagrangian describing a Fermi field coupled to a prescribed 
pseudoscalar field as a background field, in (1 + 1) dimensions. The background field is chosen to 
be in the form of the soliton of the sine-Gordon model. As is well-known, the solution of the sine- 
Gordon model is a real soliton which means that it is a solitary wave whose energy density profile is 
asymptotically in time restored to its original shape and velocity after collision with another soliton. It 
is important to note that in our model it is not the solution of E-L equation and is just a background 
field that does not have any dynamics. However, we know that the back-reaction of the fermion on 
the soliton is small and it changes the shape of the soliton slightly [211 [23 ■ Thus, it inspired us to 
investigate this coupled fermion-soliton model with sine-Gordon soliton as a background field. Since 
the equations of motion of the fermion for our model are not analytically solvable, we solve them 
using numerical methods and obtain both the bound and the continuum states for the Fermi field. Our 
purpose is to calculate the Casimir energy for this system and look closely into the details of this energy 
for our system. To this end we obtain the Casimir energy of the positive and negative parity eigenstates 
individually and then we explore these results by comparing them with the bound energy levels of the 
fermion. For the models with parity symmetry and linear equation of motion, it is possible to find 
parity eigenstates by mixing independent scattering states in a way that respects parity condition. 
Each bound energy level has specific parity. Therefore, it is worth to try to find parity eigenstates for 
continuum states to look closely the system in detail. In this paper, we compare our result with a a 
similar model in which the prescribed pseudoscalar field has a simple piece-wise form which behaves 
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like a soliton in (1 + 1) dimensions and the simple form of it makes the problem exactly solvable [28 . 
We computed the Casimir energy of this simple exactly solvable model (SESM) by subtracting directly 
the vacuum energy of the system in the presence from the absence of the disturbance |25) . However, 
it is not possible for the present model, since it is not exactly solvable and it needs another indirect 
method to calculate the Casimir energy. We take advantage of the phase shift method to calculate it. 
Moreover, the comparison between these two models helps us to investigate how the functional form 
of the background field and the details of its shape affect the result for each parity and the scattering 
ones. Meanwhile, we check the consistency of the Levinson theorem for the phase shift of each parity 
eigenstate and also the scattering states. 

In section II we introduce the model. In section III we explain how to find the bound and continuum 
states and plot some examples of them. In section IV we focus on the calculation of the Casimir energy 
by the use of phase shift method. Also, we consider the consistency of the resulting phase shifts with 
the weak and strong forms of the Levinson theorem. We then show the Casimir energy as a function 
of the parameters of the pseudoscalar field. In section V we consider a system consisting of a valence 
fermion present in the ground state, to investigate the stability of the sytem. In all the sections we 
compare the results with the ones of SESM. Finally, in section VI we conclude with the brief summary. 

2 The model of a fermion in the presence of the prescribed 
sine-Gordon 

We couple a Fermi field with a prescribed pseudoscalar field, in the form of the soliton of the sine- 
Gordon system, through the following Lagrangian 

C = $ (irfdn - Me^^ 5 ) V, (1) 

where the external field <fr(x) is chosen to be <f>{x) — (m/y/\) tan -1 [Exp(mx)] — 7rm/VA. The parame- 
ters M and m refer to the masses of the Fermi and pseudoscalar fields, respectively. Our purpose is to 
compute the Casimir energy of this system using the spectrum of the Fermi field. The Casimir energy 
in the standard approach is given by computing the change in the zero-point energy of the field due to 
the presence of the boundary conditions or the background field. This energy for a Fermi field can be 
written as follows 

/+oo p-\-oo j r-\-oo 

oo dx j o £ E (-VF+m*) 44 + J ^ d*£ (£fc und ) xl hi x2 hi 

j — i i 

/+oo r+oo ii 

/•+<» a# 

= Etd - / dk^TAP (p sea (fe) - P r(k)) + T - (2) 
Jo * 

The first equality is the relation we derived for the Casimir energy of a Fermi field in |25) . The 
functions v 3 v (x) and Vk (x) refer to the normalized fermion wavefunctions for the continuum states with 
negative energy in the presence and absence of the disturbance, respectively. The functions X2t>; (%) 
are the normalized fermion wavefunctions for the discrete bound states with negative energy and 
Abound denote the energies of these negative bound states. The extra M/2 in the last line takes into 
account the contribution from the half-bound fermionic state at E = —M in the free case. The factor 
(p sea (fc) — Po ea (k)) is the difference between the density of the continuum states with the negative 
energy in the presence and absence of the disturbance. 

If for a system we have all the normalized continuum wavefunctions with negative energy in the 
presence of the disturbance, vi{x), and energies of all the negative bound states, E^ und , we can in 
principle obtain the Casimir energy using any of the two equalities in Eq. @, though the calculations 
could be very difficult and we might need to calculate the integrals numerically. In a previous paper 



3 



[2"5] we considered a model described by the Lagrangian ([1} in which the pseudoscalar field 4>{x) is 
prescribed and in the following form 



<t>(x) 



—9q for x —I, 
fix for — Z ^ x ^ I 
+9q for Z $C x. 



(3) 



'o 





V, 1 



slope = yLf 



'o 



Figure 1: The solid and dashed lines show <j>(x) for the soliton of sine-Gordon and the simple exactly solvable 
model, respectively. The parameters #o and /j, are also shown in the figure. 

This form of 4>(x) along with the soliton of the sine-Gordon is shown in Fig. ([I]). In this figure we 
indicate the parameters 0q and /x, which are the value of 4>{x) at spatial infinities and the slope of 
the pseudoscalar field at x = 0, respectively. The form chosen for <f>(x) in Eq. ^ renders the problem 
exactly solvable and we are able to obtain the whole spectrum of the fermion as has been done in 
[281. Therefore, for this problem we could calculate the Casimir energy by the use of the first equality 
in Eq. ([2}. However, in the present problem, the form chosen for the pseudoscalar field makes the 
equations of motion analytically unsolvable. We can obtain the spectrum of this system by the use of 
an appropriate numerical method. 

It is worth noticing that the expression in the first equality of Eq. @ only contains the negative 
energy states. However, in [25] we showed that although the Lagrangian (JXJ) in (1 + 1) dimensions is 
not invariant under the charge conjugation and the particle conjugation symmetries, the total energy 
density with negative and positive energy are exactly the mirror images of each other for SESM. 
Therefore, the Casimir energy for this model can be computed only from the negative states or only 
the positive states, or the average of all of the states and the results are exactly the same in all cases. 
Note that since all the symmetries of the model with the soliton of the sign-Gordon are the same as 
SESM, the aforementioned argument is also true for the present model. 

3 The spectrum of the fermion 

The spectrum of a fermion can be totally distorted due to its coupling to a background field. In 
the presence of the background field the spectral deficiencies can be observed in the continua, and 
also the bound states can appear. To obtain the distorted fermionic states and their energies in the 
presence of the prescribed soliton of the sine-Gordon, we start with considering the Dirac equation of 
the Lagrangian ([T]). Choosing the representation 7 = 01, 7 1 = £03 and 7 s = 7°7 1 = a 2 for the Dirac 
matrices, this equation becomes 



ia\dtip — cr^dxip — M [cos <j> {%, t) + 102 sin</> (x, t)] ip = 0, 



(4) 
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where ip 



i>2 



Using the definition £(x,t) = e 



-iEt 





( i>i 


+ #2 \ 


V U*) ) 


~ U 


-#2 / 



obeyed by and £2(2;) can be written in the following matrix form 



id x — E 



—id x — E 



the equations 



(5) 



To obtain the fermionic bound states and their energies, we use a numerical method called the relaxation 
method, as we did for a similar problem in an earlier paper |27j . This method is used for solving the 
boundary value problems and it determines the solution by starting with a guess and improving it, 
iteratively. We refer you to [2 7) to become familiar with the procedure, in detail. The only difference 
is that in the present problem here the pseudoscalar <p(x) is prescribed and without any dynamics, but 
4>{x) in [57] is not prescribed and is allowed to find its final form through the equations of motion. 
Therefore, the number of the coupled first-order ODEs in the present problem is less by two. In our 
problem we use the fermionic bound states and their energies for the exactly solvable model, obtained 
in [28j . as a good initial guess. From now on we rescale all the parameters with respect to the mass of 
the fermion (M), for simplicity. 

Figure @ shows some examples of the bound state energies obtained from the numerical results. 
The left graph shows the bound energy levels of the fermion as a function of n at 9o = it, i.e. a soliton 
with winding number one and the right graph shows the bound energies as a function of 8q when ji = 1 0. 
In both graphs we also depict the bound energy levels of the fermion for SESM, for comparison. The 
solid and dashed lines refer to the energies of the model with the soliton of sine-Gordon and SESM, 
respectively. The parity of each bound state is denoted by ± signs. 





Figure 2: The energies of the bound states of the fermion. The left graph shows the bound energies as a 
function of /x at 8o = tt and the right graph shows the bound energies as a function of 8o at /i = 10. In both 
graphs the solid and dashed lines are for the model with the soliton of the sine-Gordon and SESM, respectively. 
The ± signs indicate the parity of each bound state. 



Now, we focus on the continuum states. We first obtain the fermion continuum wavefunctions for 
a scattering process in which a plane wave enters the scattering region from the left. Then, using 
the scattering wavefunctions, we find the parity eigenfunctions which satisfy the parity condition 

Pt(x,t) = -V 2 {(x,t). 

Choosing £i(x) = r)i(x) + ir]2(x) and ^(x) — f}a{x) + iri4,{x) in which in general rjiS are real functions 
of x, we obtain the equations of motion satisfied by rjiS in the following form 



r] 1 + cos <j) {x) f]3 — Er]2 - sin <j> (x) 774 = 0, 
r\ 2 + cos <f) (x) 174 + Erfi + sin <j> (x) i] 3 = 0, 
?7 3 + cos (j) (x) r)i + E7J4 + sin <j> (x) i] 2 = 0, 
f] A + cos (f> (x) r/ 2 - Er) 3 - sin (x) rji = 0. 



(0) 
(7) 
(8) 
(9) 
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For the form chosen for <f>(x), i.e. the soliton of the sine-Gordon, we cannot solve this set of equations, 
analytically. Therefore, we use an appropriate numerical method to find the scattering wavefunc- 
tions. We solve this set as an initial value problem and find the wavefunctions for the whole interval 
(— oo,+oo), using the so-called Runge-Kutta methods. In order to find the solutions with high accu- 
racy, we use the Runge-Kutta-Fehlberg method of order 6. Since the form chosen for 4>{x) in SESM is 
similar to the soliton of the sine-Gordon as x — > ±00, we can take advantage of the solutions of SESM 
to determine the initial boundary values for solving the equations. We already have all the solutions 
of this model, including the scattering wavefunctions [28] Figure © shows r)iS for the case with 
the parameters 60 = 7T, h = 10, k = 3.0 and E = +y/k 2 + M 2 . This figure also shows r^s of the SESM, 
for comparison. 
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Since the set of equations ([5]) is a linear set for the Fermi field, a linear combination of its solutions 
is another solution. We take advantage of this fact to find the continuum parity eigenf unctions. We 
choose two independent scattering solutions <pk(x) and 4>-k{x) obtained by the numerical calculations 
mentioned above. We linearly combine these two solutions and determine the coefficients such that the 
combined solution becomes the simultaneous eigenstate of the Hamiltonian and the parity operators. 



In Figs. Q and ([5]) we depict the parity eigenstates for the case with the parameters 9q — w, [i = 10, 
k = 3.0 and E = +yjk 2 + M 2 . In this figure we also show the continuum parity eigenstates of 

SESM with dashed lines. ( V \ [ X \ + %V \ [ X \ ] and ( V l j X ) + W2 _ [ X \ ) arc the continuum parity 

eigenstates with the positive and negative parity, respectively. The relations y^(—x) = —y^(x), 
vt(~ x ) = ut { x )t Hi (~ x ) — Vi{ x ) an< i V2{~ x ) = — y-i( x )i which are equivalent to the parity operation 
P£(x,t) = —(72^(x,t), can be easily checked in these graphs. 
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4 The calculation of the Casimir energy using the phase shift 
method 

Now, having all the spectrum of the system, we can calculate the Casimir energy and explore it in 
detail. The continuum parity eigenstates obtained for our problem are not normalized. Therefore, we 
cannot directly use the relations given in Eq. ^ to compute the Casimir energy. However, we can 
utilize the phase shifts of the continuum states to calculate this energy. 

4.1 The phase shift and Levinson theorem 

Our starting point is the following relation 

p(k)-p (k) = -^-6(k), (10) 

where p(k) — po(k) is the difference between the density of the continuum states in the free and the 
interacting cases and 8(k) = S sky (k) + S sca (k), i.e. 5(k) sums over the contributions from both positive 
and negative energies. This relation also holds for the sea and sky, separately. 

Substituting Eq. (fTU)) into the second term in the relation of the Casimir energy, Eq. and then 
integrating by parts, we obtain 

+oo /-+OC 11 J 

dkVk 2 + M 2 {p sea (k) - /9 s oa ( fc )) = - / — Vk 2 + M 2 — (<5 sea (fc) - 5 sca (oo)) 

Jo tt dk 

j=== {5 sca (k) - <5 sca (oo)) + —M (<5 sca (0) - <f ca (co)) . (11) 

In the first equality we have just subtracted a zero term (^<5 sca (+oo)) from the original one. Now, we 
can compute the second term in the expression of the Casimir energy using the phase shift. 

Comparing the coefficients of e lkx on the left- and right-hand sides of the scattering region, we could 
obtain the scattering matrix element, which is related to the scattering phase shift as S(k) = e 15 ^. As 
one knows, the parity eigenfunction in both limits x — > ±oo is a linear combination of e lkx and e~ zkx . 
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Figure 5: The wavefunction of the fermionic state with negative parity. The graphs of y~ (x)s as functions of 
the spatial variable x, for the parameters do — it, fi = 10, k — 3.0 and E — +\/k 2 + M 2 . The solid and dashed 
lines show the normalized rji(x)s for our model and SESM, respectively. 



Thus, the scattering matrix clement for parity eigenstates is obtained by comparing the coefficients of 
e %kx Q f |_ ne p ar ity eigenfunctions on the left- and right-hand sides and is related to the phase shift as 
S±(k) = e M ±W. 

Now we check the consistency of the resulting phase shifts with the Levinson theorem. The weak 
form of this theorem for the Dirac equation is as follows |H] 

AS = [<f ky (0) - <f ky (oo)] + [<P oa (0) - <P oa (oo)] = ( N + -1 - -± J 7T, (12) 

where N is the total number of bound states, including positive and negative ones, N t the total number 
of the threshold bound states at the given strength of the potential, and N® the number of bound states 
at zero strength of the potential, i.e. the free Dirac case. The strong form of the Levinson theorem 
relates the value of the phase shift at each boundary of the continua to the number of levels which have 
crossed those boundaries, in the process of building up the disturbance. This form of the theorem can 
be expressed in the following form for k = 

S{0) = (N cxit - AUer) 7T. (13) 

This relation holds for each of the continua, separately. For each continuum iV ex it (N ontcr ) is the number 
of the bound states that exit (enter) that continuum from that boundary (E — +1.0 or E — —1.0) 
as the strength of the potential is increased from zero to its final finite value. In this equation the 
threshold bound states, mentioned above, should be included in A^it and N cntcr , with the coefficient 
1/2. Moreover, the strong form of the Levinson theorem for k = oo can be written in the following 
form 

5(oo) = (N enteI - iV exit ) 7T. (14) 

This relation also holds for each continuum, separately. For each continuum ./V cx it (N cntcr ) is the 
number of the bound states that exit (enter) that continuum from k — oo (E = +oo or E = — oo) as 
the strength of the potential is increased from zero to its final finite value. Analogous relations hold 
between each parity phase shift and the corresponding bound states as follows 

/ jy N° \ 

AS ± = [6t y (0) - St y (^)} + [*T (0) " <5± a M] = [N± + ~ U (15) 
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<5± (0) = (N cxit ,± - iV cntol%± ) 7r and <5±(oo) = (iV e nter,=t - N cxit: ±) n. 



(16) 



As an example of our results, in Fig. © we plot the phase shift for our system as a function of 
k, for the parameters 9q = ir and \i = 10. The left graph shows the 5 sky (k), i.e. the phase shift 
for the scattering states with the positive energy +\ , k 2 + M 2 and the right graph shows S^ y (k), i.e. 
the phase shift for both parity eigenstates with the positive energy. In both graphs we also depict the 
phase shift of SESM with the same parameters, for comparison. It is easy to check that the phase shifts 
depicted in these figures are consistent with both the weak and strong forms of the Levinson theorem. 
In particular i5 sky (fc — > +oo) = +(9o and 5 s ± y (k —> +oo) = +6q/2, which are consistent with the results 
of the adiabatic method of Goldstone and Wilczek [30]. For these parameters we can see from Fig. J2j 
that two levels (the first one with positive parity and the second one with negative parity) have exited 
from the sky, and 7V t °'£ ky = 1 and iV t l ky = at E = +1.0. Therefore, using the strong form of the 
Levinson theorem, we expect 5 sky (0) / 7r = 3/2, <fL y (0)/7T = 1/2 and <J!_ y (0)/7r = 1. Moreover, the 
equality S sky (k) = 5 s + y (k) + S s ^ y (k) holds for every k. 




Figure 6: The graphical representation of phase shifts of the continuum states with positive energy, for /i = 10 
and #o = 1"- The left graph: the phase shift of the scattering state, <5 sky (fc). The right graph: the phase shifts 
of the ± parity eigenstates, S J y (k) . In both graphs the results for the models with the soliton of sine-Gordon 
and SESM are shown by the solid and dashed lines, respectively. The ± signs in the right graph show the 
parity of the corresponding eigenstates. 



4.2 The Casimir energy 

Now, using the last expressions in Eqs. and (jlll) and the bound states shown in Fig. @ we calculate 
the Casimir energy for different values of the parameters of the prescribed soliton. 

In Fig. ([7]) we plot the Casimir energy as a function of fj, at 9q = i-e. a soliton with winding 
number 1, for our model and for SESM by the solid and dashed lines, respectively. The left and right 
graphs in this figure, respectively, show the total Casimir energy and the Casimir energy for the each 
parities. The parity of states in the right graph is denoted by the ± signs. As can be seen from the left 
graph, in both models, there is a sharp maximum occurring when the bound energy level crosses the 
line of E — 0. This crossing happens for a larger value of fi for our model, based on the bound energy 
levels shown in the left graph of Fig. @. For our model the bound energy level crosses the line of 
E = at /i 4.000, while for SESM this crossing occurs at the lower value of the slope, i.e. /i ~ 2.957. 
Also, the value of the Casimir energy is lower in the case of the soliton of the sine-Gordon. The largest 
difference between the graphs of these two models happens around the maximum, as is shown in the 
zoomed box of the left graph. In both models the Casimir energy reaches the expectable value of zero 
when the slope of 4>(x) at x — decreases to zero, i.e. when the vacuum energy approaches that of 
the trivial vacuum, despite the residual non-trivial boundary conditions. Also, the Casimir energy for 
both models have the same limit when the slope of the pseudoscalar field tends to infinity. This limit 
is zero at 8q = nn, i.e. when we have a proper soliton with winding number n. However, for other 
values of 9q the Casimir energy is in general non-zero, when /i — ¥ oo in both models (see Eq. (3.9) in 
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[25.). From the graphs of the energy for the separated parity states we can see that there is a sharp 
maximum for the Casimir energy of the positive parity occurring when the first energy level in the 
left graph of Fig. © with positive parity crosses the line of E = 0. However, the Casimir energy of 
the negative parity does not have any maximum, since none of the negative parity levels crosses the 
line of zero energy at 8q — ^ ■ The total Casimir energy in the left graph is the sum of the Casimir 
energy for the positive and negative parity states shown in the right graph. One can indeed calculate 
the total Casimir energy using the scattering phase shift and the whole bound states of the fermion, or 
equivalently by adding the Casimir energy obtained for each parity using the phase shift of the parity 
eigenstates and the corresponding bound states. 




0. 50. 100. 150. 200. 250. ' . 50 , | 00 . 150. 200. 250. 



Figure 7: The graphical representation of the Casimir energy as a function of fi, the scale of variation of the 
soliton, at Go — tv. The left graph shows the total Casimir energy for the model with soliton and SESM, with 
the solid and dashed lines, respectively. In the box we focus on the small values of /i to show the details of the 
maximum and the differences between the results of two models. The right graph shows the Casimir energy 
for the ± parities, separately. 



In Fig. (]8]) we present the Casimir energy as a function of 8q for /i = 10, by the solid and dashed 
lines for the model with the soliton of sine-Gordon and SESM, respectively. As before, the left and 
right graphs in this figure show the total Casimir energy and the Casimir energy for the each parities, 
respectively. As can be seen, the total Casimir energy is, on the average, an increasing function of 9q 
for both models and there are two mild cusps in each graph. Comparing these graphs with the right 
graph of Fig. ([2]), we conclude that these cusps occur when the bound energy levels cross the line of 
E = 0. The bound energy levels cross E — at 0$ w 0.6257T and 9 w 1.8757T for our model, and at 
#o ~ 0.576-7T and 9 1.5967T for SESM, in the interval ^ 6q ^ 2-7T. From the right graph it can be 
seen that the Casimir energy for both parities has a mild cusp occurring where the levels with positive 
and negative parity cross the line of E = 0, as expected. As before, the total Casimir energy in the 
left graph is the sum of the Casimir energy for ± parities in the right graph. 



5 Stability of the solutions 

In this section we consider a system consisting of a valence fermion in the ground state and explore 
the effect of the Casimir energy on its total energy. The total energy for such a system is the sum of 
the Casimir energy and the energy of the valence fermion. Note that the energy of the valence fermion 
should not be added when this energy is negative, since it has already been taken into account in the 
Casimir energy. We show the total energy in Figs. ([§]) and (fTTJl) . In the left graph of Fig. © we vary 
the slope /x of the prescribed pseudoscalar field while its asymptotic value is fixed at 9q = n and in the 
right graph of Fig. (fTTJ)) we plot the total energy for n = 10 and different values of 0q. In both graphs 
the solid and dashed lines refer to the results for the model with the soliton and SESM, respectively. 
Using the Casimir energy for the individual parities shown in the right graphs, we are able to explore 
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Figure 8: The graphical representation of the Casimir energy as a function of 9 , the value of the soliton at 
spatial infinity, at /i = 10. The left graph shows the total Casimir energy for the model with soliton and SESM, 
with the solid and dashed lines, respectively. The right graph shows the Casimir energy for the ± parities, 
separately. 



the stability of each parity, separately. It should be mentioned that since the ground state of the 
fermion has positive parity, the total energy for the negative parity has only one contribution coming 
from the Casimir energy for the negative parity eigenstates. We show the total energy for the positive 
and negative parities with the ± signs. In the right graph we can see an interesting result. For both 
models there is a minimum occurring at 8q w 7t, which corresponds to a soliton with winding number 
one. This means that not only this configuration is energetically favorable, but also it is stable against 
small fluctuations in the parameters of the background field when this field is a soliton with a proper 
winding number, as expected. 




Figure 9: The left graph shows the total energy (the sum of the energy of a valence fermion and the Casimir 
energy) as a function of /i when #o = it. The solid and dashed lines represent the total energy for the model 
with the sine-Gordon soliton and SESM, respectively. The right graph shows the total energy for positive and 
negative parities, separately. 



6 Conclusion 

We compute the Casimir energy for a coupled fermion-pseudoscalar field in which the pseudoscalar field 
is prescribed and chosen to be the sine-Gordon soliton. This model is not exactly solvable. Therefore, 
we use the relaxation method and the Runge-kutta-Fehlberg method to obtain the bound states and 
continuum scattering wavefunctions of the fermion, respectively. Since the model is invariant under the 
parity and its equations of motion are linear, the parity eigenfunctions of the system can be obtained 
using a linear combination of the continuum scattering states. The coefficients of this combination 
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Figure 10: The left graph shows the total energy (the sum of the energy of a valence fermion and the Casimir 
energy) as a function of 80 when fi = 10. The solid and dashed lines represent the total energy for the model 
with the sine-Gordon soliton and the SESM, respectively. The right graph shows the total energy for positive 
and negative parities, separately. Notice that the system attains its lowest energy at 80 = 7T. 

arc determined by the use of the parity condition. It is interesting to find the continuum parity 
eigenfunctions, because bound states are eigenstates of the parity. Having all these states, it is possible 
to investigate the system and all the results for each parity, in detail. To end up with Casimir energy, 
we use the phase shift method. To be sure about the resulting phase shifts, we check their consistency 
with the weak and strong Levinson theorems for each parity and scattering phase shifts. Moreover, we 
compute and plot the Casimir energy as a function of the slope of the pseudoscalar field at x = (/z) 
and its value at infinity (9o). The graph of the Casimir energy as a function of \i has a sharp maximum 
at the value of fx in which the fermion has a zero-mode and goes to zero for /i — > 0, and also for /i — > 00 
for a soliton with the proper winding number. In the graph of the Casimir energy as a function of 6*o, 
there are some cusps when the fermion bound energies cross the line of E = 0. Also, we calculate and 
plot the graph of the Casimir energy as functions of /i and 6$ for each parity to investigate what exactly 
happens in the system by comparing them with the fermionic bound energy levels. The Casimir energy 
obtained by the use of the scattering phase shifts is the sum of the Casimir energy obtained from both 
parity phase shifts. Finally, we study the stability of the system in the presence of a real fermion in 
the ground state. We show that the total energy has a minimum, independent of the details of the 
shape of the background field, for a proper soliton with winding number 1. For each of the graphs we 
compare the results with SESM that mimics the soliton in (1 + 1) dimensions. 
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